Unusual Double-peaked Emission in the SDSS Quasar 
J093201.60+031858.7 



R.S. Barrows a '*, C.H.S. Lacy a < b , D. Kennefick a < b , J. Kennefick a ' b , M.S. Seigar c 

a Arkansas Center for Space and Planetary Sciences, 202 Old Museum Building, University of Arkansas, 

Fayetteville, AR 72701 

b Physics Department, University of Arkansas, 825 West Dickson Street, Fayetteville, AR 72701 

c Department of Physics and Astronomy, University of Arkansas at Little Rock, 2801 S. University Avenue, Little 

Rock, AR 72204 



Abstract 

We examine spectral properties of the SDSS quasar J093201. 60+031858. 7, in particular the presence 
of strong blue peaks in the Balmer emission lines offset from the narrow lines by approximately 
4200 km s _1 . Asymmetry in the broad central component of the H/3 line indicates the presence of 
a double-peaked emitter. However, the strength and sharpness of the blue H/3 and blue H7 peaks 
make this quasar spectrum unique amongst double-peaked emitters identified from SDSS spectra. 
We fit a disk model to the H/3 line and compare this object with other unusual double-peaked quasar 
spectra, particularly candidate binary supermassive black holes (SMBHs) . Under the binary SMBH 
scenario, we test the applicability of a model in which a second SMBH may produce the strong 
blue peak in the Balmer lines of a double-peaked emitter. If there were only one SMBH, a circular, 
Keplerian disk model fit would be insufficient, indicating some sort of asymmetry is required to 
produce the strength of the blue peak. In either case, understanding the nature of the complex line 
emission in this object will aid in further discrimination between a single SMBH with a complex 
accretion disk and the actual case of a binary SMBH. 
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1. Introduction 



Quasars and active galactic nuclei (AGN) with spectra showing double-peaked emission lines are 
an interesting class of objects, though in many cases the origin of the double peaks is not fully 
understood. Some of these spectral features have been examined in the context of emission from an 
accretion disk. These objects are commonly kn own as double-peaked emitters (or di s k emitters) and 



have been the subjects of a number of studies (jEracleous and Halpernl . 



1994, 



2003; 



Strateva et al 



2003). Double-peaked emitters are important objects to study because the double-peaked profiles 



provide information about the structure of the accretion disk around supermassive black holes 
(SMBHs). The disk profiles are typically centered around Balmer emission line p ositions presumably 
becau se of photoionization and electron scattering in a geometrically thin disk (jChen and Halpernl . 



1989). In the disk profiles, there is typically a peak on the red side and a peak on the blue side of 



the central emission line, though in some cases the peaks are very broad and shallow. Furthermore, 
the blue peaks are typically stronger than the red peaks because of relativistic beaming. 

For many double-peaked emitters, a simple (circular, Keplerian) disk model provides a satisfactory 
model. However, this sub-class of AGN contains a range of Balmer emission line profiles, and in some 
cases the fits require additional complexities (i.e. asymmetries) such as eccentricity in the disk or a 
hot spot on the disk surface. In a few cases, only one peak appears to be offset from the narrow lines, 
and such objects have been interpreted as candidate binary or recoiling SMBHs because the spectra 
resemble multiple emission line regions. These interpretations are important to consider because 



binary and /or recoiling SMBHs are expected to result from galaxy mergers (jMerritt et al 



Komossa , 



2004; 



2006) and therefore many such cases should exist. Additionally, they are believed to play a 



role in the evolution of the SMBH population ([Hopkins et al 



20051 ). Therefore, there has been much 



interest in the possibility of identifyi ng close binary SMBHs by locatin g double-peaked emission lines 



in the spectra of AGN and quasars (Gaskcll, 



Gaskell . 


1996; 


Zhou et al.. 


2004) 



20041 ). In these cases, determining 



the true physical origin of the double-peaks is difficult since they may be partly or entirely due 
to complex disk emission instead. Interestingly, a few candidate binaries have been identified. 
For example, the object SDSS J153636. 22+044127 . (fro m here on SDS SJ1536+0441) has bee n 
interpreted as a binary SMBH (Boroso n and Lauerl 12009). a disk emitter (jChornock et all 120101 ) . 



and as both (jTang and Grindlav . 



20091 ) . 



Another candidate, SDSS J10504 1.35+345631.3 (from 



here on SDSSJ1050+3456) shows similar features and was first identified by 



Shields et al 



(120091 ) 
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as a possible recoiling or binary SMBH. A third candidate, SDSS J092712. 65+2943 44.0 (from here 



on SD SSJ0927+2943), was first identified and inter preted as a recoil i ng SMB H (Kom ossa et al 



20081 ). and more recently also as a binary SMBH (|Bogdanovic et al 



2009b 



Dotti et al 



20091 



The diversity of these interpretations indicates that accretion processes around a singe SMBH 
or in a binary SMBH are expected to produce a range of complex spectral signatures. Results 
from sim ulations indicate that profiles of close binaries strongly resemble those of double-peaked 



emitters (jBogdanovic et al 



2008 



2009al) ; therefore, it is important to examine double-peaked 



profiles, particularly those with strong and distinct multiple peaks, as evidence for complex accretion 
processes that may be the result of a binary SMBH. 

Inspired by the above mentioned previous searches, and the results of those searches, we performed a 
similar search through the Sloan Digital Sky Survey (SDSS) archives to identify spectra with double- 
peaked emission lines at z < 0.89; this redshift limit was chosen so that all spectra would include 
the H/3 line. All matches were individually inspected. In this paper, we describe the spectrum of a 
quasar, SDSS J093201.60+031858 (from here on SDSSJ0932+0318), that has a blue set of Balmer 
emission lines offset from the narrow lines by approximately 4200 km s" 1 . While the H/3 line has a 
broad, red peak that is indicative of a double-peaked emitter, it has unusually strong blue H/3 and 
H7 peaks that are unusual among SDSS identified double-peaked emitters. We examine similarities 
and differences between the features of SDSSJ0932+0318 and those of other unusual double-peaked 
emitters, particularly candidate binary SMBHs, to help determine the physical nature of the line 
emission. Throughout the paper we adopt the cosmological parameters H Q = 71 km s _1 Mpc -1 , 
VLm — 0.27 and Cl vac — 0.73. 



2. Spectral Analysis 



2.1. The Continuum 



The quasar spectrum w e discuss was drawn from the fourth edition of the SDSS quasar catalogue 



(Schneider et al 



20071 ). It has SDSS point source magnitudes of u = 18.731, q = 18 460 



18.271 and z = 17.945. Using the suggested conversion from 



Jester et al 



18.386, 



(120051) . B = g 



0.17(w — <?)+0.11, yields B — 18.39, which is also corrected for a Galactic colo r excess of E(B — V) 



0.053, taken from the NASA Extragal actic database follow ing the model of 



Schlegel et al 



Converting to a bolometric luminosity ([Marconi et all 120041 ) and correcting for a stellar contribution 



(1998). 
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of 33 percent (the mean value from lEracleous and Halpernl . 120031) yields Lbol = 4.25 x 10 erg 



s _1 . The spectrum was corrected for Galactic extinction using the IRAF task 'deredden.' Since 
disk emission is evident in the spectrum (figure 1), the line of sight to the central SMBH is probably 
not significantly obscured so we did not correct for internal extinction. The continuum slope was 
modeled as a power law (fx = A~ Q ), and the best fit continuum slope has an index of a = 2.05. 

2.2. The Disk Profile 

The Ha line has been redshifted almost entirely out of the spectrum, therefore we primarily focus 
on the H/3/[OIII] complex. To fit the asymmetric line profile and to test the applicability of the 
disk-emitter scenario, we fit a disk profile to this region. The disk model was centered at the redshift 
of the [OIII]5007 line based on the assumption that it represents the re dshift of the central SMBH 



(discussed in section 3.2). The disk model we used was adopted from IChen and Halpernl (J1989) 



and assumes a circular, Keplerian disk. The assumption of circular shape is used since an accretion 
disk, in the absence of external influences, should be circularized due to the energy loss within the 
disk. In the model, which has been widely used to model double-peaked profiles, the parameters 
are as follows: the local broadening of the line (a), the inclination of the disk (i), the inner and 
outer radii of the disk (£i and £2) in gravitational units {r g ), and emissivity as a function of radius 
on the disk (e = where q is a dimensionless power law index. In previous uses of this model, 
the power law index q is usually constrained near 3. Therefore, we have set q = 3 to minimize the 
number of free parameters. The local broadening tends to range from a ~ 500 — 2000 km s -1 , and 
we used a value of 1000 km s _1 which is within that range and typical of disk model fits in the 
literature. The inclination, inner radius and outer radius are therefore the three parameters that 
we adjusted to find the best fit, and each best fit was obtaine d by eye with the restricti on that the 



fit be bounded by the observed emission line as was done in IChen and Halpernl ([19891 ) . We fixed 
the inclination at three values, 15°, 30°, and 45°, and varied the inner and outer radii to obtain the 
best fit for each inclination. For 15°, the best fit produces a sharp red peak that does not match 
the broad, shallow red peak that is observed. In contrast, for 45°, the best fit does not produce 
the extended red wing that is blended with the [OIII] emission lines. We adopted an inclination of 
i=30° and inner and outer radii of £1 = 225 r g and £2 = 2700 r g , respectively, for which the broad 
and shallow red peak is best fit (figure 1). This is important since, after subtraction of the disk 
and subsequent fitting of the emission lines (section 2.3), the [OIII] 5007/ [OIII] 4959 ratio is ~3, in 
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Figure 1: H/3/[OIII] complex at the rest wavelength with the disk profile overlaid. The central wavelength of the 
disk model is set to the position of the narrow H/3 line. The parameters of the fit are as follows: q = 3, a = 1000 km 
s -1 , i = 30°, §i = 225 r g and §2 = 2700 r g . Notice the strong blue peak that is not adequately fit by the circular, 
Keplerian disk model. 
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agreement with atomic theory. The correct position of the blue peak is also produced when these 
parameters are used. However, in none of the fits does the symmetric disk profile account for all 
of the flux in the blue peak which is unusually strong for a double-peaked emitter. A similar fit 
was performed on the H7 region. Though the H7 emission is weaker and a fit is more difficult, it 
is clear that a symmetric disk profile is not sufficient there for the same reason that the blue H/3 
peak is not well modeled. 

2.3. Emission Lines 

After subtraction of the disk profile, Gaussian profiles were fit to the broad and narrow central 
components of the H/3/[OIII] and H7 regions using the IRAF task 'fitprofs.' In the H/5/[OIII] 
region all parameters of the fit (peak position, peak flux and FWHM) are left as free parameters, 
but the peak position and FWHM in the H7 region are fixed because the features are much weaker. 
The FWHM of the narrow lines ([OIII]5007, [OIII]4959, H/3, H 7 , [NeIII]3869 and [OII]3727) range 
from 300 - 700 km s -1 and probably originated in the same narrow line region (NLR). The FWHM 
of all the broad central Balmer components (H/3 and H7) are also similar, indicating that they 
originated in a region nearer the central black hole, typical of a broad line region (BLR). Gaussian 
profiles were also fit to the residuals of the blue H/3 and H7 peaks. Table 1 lists the redshifts, 
Gaussian FWHMs and the line flux ratios relative to H/3 for all of the fits. 



3. Interpretation 

3.1. Comparisons with Candidate Binary SMBHs 

While the asymmetry of the broad H/3 component of SDSSJ0932+0318 that is blended with the 
[OIII] lines is likely the red shoulder of a broad double-peaked emitter, the strength of the blue 
H/3 p eak relative to the red peak is unusual for double-peaked emitters ([Eracleous and Halpernl . 
19941 ). Additionally, the blue peak is offset from the central narrow line by ^4200 km s _1 which 
is on the high end fo r double-peaked emit ters (larger than 91 percent of the sample of 116 double- 



peaked emitters from 



Strateva et al 



2003) . These properties warrant comparison with other objects 



identified from SDSS spectra that exhibit unusual double-peaked line profiles and which have been 
proposed as candidate binary SMBHs (figure 2). 



G 



SDSSJ1536+0U1: The extended red wing and strong blue peak of SDSSJ0932+0318 is similar 
to the features seen in SDSSJ1536+0441. The blue peak in SDSSJ1536+0441 is blueshifted by 
^3500 km s , and the unusually high ratio of the blue peak to the broad, red peak makes it an 
extreme double-peaked emitter. As such, SDS SJ1536+0441 has re c ently been proposed as both 



a double-peaked emitter and a binary SMBH ([Tang and Grindlav 



2009). This last hypothesis 



is particularly intriguing since it suggests that strong peaks in some double-peaked emitters may 
actually be caused by a second, smaller SMBH. In this scenario, the less massive secondary SMBH 
is orbiting the primary SMBH and moving toward us, and the broad blue Balmer components 
in the spectrum may come from an accretion flow onto the secondary, or alternatively from the 
portions of the primary a nd circum-binary accretion disks near to and heated by the secondary 



(|Tang and Grindlav 



2009) 



SDSSJ1050+3456: A similarly extreme difference between the blue and red peaks is seen in 
SDSSJ1050+3456, and could likewise be interpreted as evidence for a double-peaked emitter in 
a binary SMBH. Interestingly, however, the extended red wing is only apparent in the Ha line of 
SDSSJ1050+3456, while the broad H/3 peak appears to be symmetric. 

SDSSJ0927+2943: The SDSS spectrum of SDSSJ0927+2943 is unique in that both the blue and 
red emission line systems contain narrow line components, including in the high ionization forbidden 
lines; and the broad components, which are only seen in the Balmer and Mgll lines of the blue 
system, all appear symmetric with no visible disk profile. If interpreted as a binary SMBH, the 
red system is associated with the NLR that surrounds both SMBHs and is at the redshift of the 
primary, more massive SMBH, and the blue system is associated with the secondary, smaller SMBH 
that provides illumination for both emission line systems. 

3.2. The Binary SMBH Model 

Testing these models on SDSSJ0932+0318 and other strongly double-peaked emission line pro- 
files, is particularly important to do since spectra with double-peaked emission lines is a potential 
method of identifying binary SMBH candidates. That the profile of a circular disk is inadequate 
suggests the possib i lity o f an external influence, perhaps from a secondary SMBH. Arguments in 



Tang and Grindlav! (|2009h (and references therein) suggest that the presence of a secondary SMBH 



is effective at producing the strong peaks seen in some double-peaked emitters (and making a disk 
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Figure 2: Top left: SDSS spectrum of SDSSJ0932+0318 with double-peaked H/3 and H7 emission lines. Top right: 
SDSS spectrum of SDSSJ0927+2943 with double-peaked H/3, H7 and narrow emission lines. Bottom right: SDSS 
spectrum of SDSSJ1050+3456 with double-peaked Ho, H/3 and H7 emission lines. Bottom right: SDSS spectrum 
of SDSSJ1536+0441 with double-peaked Ho, H/3 and H7 emission lines (the right side of Ho, which is cutoff in 
the SDSS spectrum, reveals a broad, shallow red wing IChornock et all I2010P) . All the spectra have been shifted 
to the rest-frame. Notice that the double-peaks in the Balmer emission lines have varying degrees of strength and 
asymmetry. For clarity, the spectra are smoothed to 600 km s —1 . 
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profile more visible by enhancing the illumination of the disk surface), especially those that re- 
quire asymmetric disk profiles, and that many such double-peaked emitters may in fact be binary 
SMBHs. 

The models of potential interactions between binary SMBHs, accretion disks and circumbinary 
disks for the previously discovered candidates suggests that signatures of binary SMBHs are likely 
to be varied. The visibility of a disk profile in the low-ionization lines is thought to be dependent 
on a number of factors: illumination (internal and external) of the geometrically thin portion of the 
disk, the disk radial size, and the orienta tion of the disk pl a ne to t he li ne of sight. As for a second 

(120091 ) suggest 



Dotti et al 



set of high ionization lines, arguments in iBogdanovic et al.l ((2009b) and 
that the secondary SMBH may be able to produce its own NLR by ionizing the low density region 
within the circumbinary disk. The subject of our analysis, SDSSJ0932+0318, is clearly within the 
class of objects which show an apparent disk profile and for which there is no second set of high 
ionization lines. In the binary SMBH scenario, this most likely implies that the secondary SMBH 
is not powerful enough to produce emission lines with high ionization potentials and low critical 
densities such as [OIII]5007, [OIII]4959, [NeIII]3869, [OII]3727 and [NeV]3426. The observed red 
NLR is the rarified medium that is large enough to surround two SMBHs that are separated by 
one parsec or less, and it is clearly at the redshift of the primary SMBH since it also coincid es with 
the central wavelength of the disk model. Following the model of iTang and Grindlavl (|2009l ). if the 



secondary SMBH is massive enough, then it may be able to maintain and illuminate its own BLR. 
Alternatively, the seco ndary SMBH may t idally distort, or heat the outer edge of, the primary 



20091 ) . Both scenarios can explain the observed blue Balmer 



SMBH's accretion disk (ILodato et al 
peaks (figure 3). 

Binary Parameters: 

Under the binary SMBH hypothesis, it is instructive to estimate the physical parameters of the 
system to determine if such a system may be identified from multiple observations and in what stage 
of binary SMBH evolution it would be placed. The individual masses may be estimated if we make 
the assumption that the central and blue broad H/3 components arc from BLRs around the primary 
(Mi) and secondary (Mia ) SMBHs, respectively. In this case, FWHM#^ and L5100 may then be 



used with the relation of 



Vestergaard and Peterson! (|2006l ). To determine the optical luminosities, 



the total bolometric luminosity calculated in section 2.1 (including the correction for 33 percent 
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Figure 3: H/3/[OIII] complex at the rest wavelength with the disk profile plus best fit Gaussian components over- 
laid. The Gaussian comp onents were fit after subtraction of the disk model. This analysis is similar to that of 
iTang and Grindlavl i2009l V 
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Emission Line 


Redshift 


FWHM (km s- 1 ) 


Flux relative to Hf3 


Central Narrow Components: 


[OIII]5007 


0.4199 


560 


5.5 


[OIII]4959 


0.4198 


660 


1.9 


H/3 


0.4201 


310 


1.0 


H 7 


0.4197 


320* 


0.5 


[NeIII]3869 


0.4190 


1800 


2.4 


[OII]3727 


0.4197 


430 


2.5 


Central Broad Components: 


H/3 


0.4177 


2710 


1.0 


H 7 


0.4179 


2710* 


0.3 


Blue Components: 


H/3 


0.4001 


2050 


1.0 


H 7 


0.3999 


2050* 


0.4 



Table 1: Redshifts, FWHMs and line ratios for the Gaussian fits to the spectrum. The redshifts were calculated 
using the SDSS vacuum wavelengths. The line widths were corrected for an instrumental resolution of 150 km s . 
* Fixed values. 



stella r contribution) was converted to a 5100A luminosity using L5100 ~ Lbol/9 ( Peterson et al 



2004). Then, the relative fraction contributed by each SMBH was calculated using the ratio of the 
H/3 luminositi es from the line fitting and converting to a ratio of the 5100A luminosities using the 
relation from iGreene and Hoi (|2005l) . The resulting luminosities are L5100 = 2.92 x 10 44 erg s _1 
and L5100 = 2.39 x 10 43 erg s _1 for the primary and secondary SMBHs, respectively, providing 
mass estimates of Mi m 3 x 10 7 M Q and M2 ~ 5 x 10 6 M Q and a mass ratio of M2/M1 w 
0.16. Additionally, using the combined mass, and assuming they are in a circular orbit with a 
relative velocity difference of 3700 km s _1 (the separation between the central and blue H/3 broad 
components) the calc ulated separation is ^0.01 pc or ~14 1-d. The velocity difference is very similar 



to that calculated bv lBoroson and Lauerl (|2009| ). but the separation is nearly an order of magnitude 
smaller because of the smaller total mass. Though these are only rough estimates, they show that, 
if the system is a binary, it is in the final stage of the binary's lifetime. 

If the separation is small enough, there will likely be complex interactions between the secondary 
SMBH and the accretion disk of the primary SMBH. We may estimate the physical disk size using 
the primary SMBH mass, which yields inner and outer radii of ~0.0003 and ~0.004 pc, respectively. 
This result is interesting since it implies, based on the above calculations, that the secondary SMBH 
would be near the edge of the primary disk, potentially close enough to heat and /or distort it. For 



the primary SMBH, we calculated the radius of the BLR (from the relation of lKaspi et al 



20051 ) to 
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be ~0.008 pc or ~9.6 1-d. This result further implies that the binary system would be nearly inside of 
the primary's BLR. Ho wever, the scaling relations u sing AGN luminosity may not directly apply to 
double-peaked emitters ( Lewis and Eracleous , 20061 ) so these masses and orbital parameters are not 
necessarily realistic; they merely indicate that the sizes of accretion disks, BLRs and binary SMBH 
separations are of similar scales, allowing for complex interactions that are capable of producing 
the observed spectrum of many unusual double-peaked emitters. For instance, the central broad 
components are slightly blueshifted (^500 km s _1 ) relative to the central narrow components, 
indicating that there may be some sort of interaction between the secondary SMBH and the BLR 
of the primary SMBH. In this case, the broad components will not be entirely due to virialized gas 
around the SMBHs. 



3.3. Chance Projection 

Since the analysis of this quasar spectrum considers the possibility of a binary SMBH, we also 
consider the possibility of a chance projection of two quasars. However, there is a small statistical 
likelihood of a positional coincidence. Furthermore, there are no blue components to the strong 
narrow forbidden lines ([OIII] and [Oil]) that are typically seen in quasar spectra. So, if the blue 
components in the spectrum are the result of a second SMBH, the second SMBH does not have a 
typical NLR and is more likely to be within the same NLR as the primary SMBH. 



3-4- Alternative Mechanisms for Asymmetry in the Disk 

In this section we investigate the possibility that only one SMBH is present and that the double 
peaks are produced entirely by a single accretion disk. Based on the disk model fit in figure 
1, this requires some form of asymmetry in the disk. The structure of accretion disks around 
SMBHs has b een widely studied, with various models proposed. The model for the disk fitting 
we used (from I Chen and Halpernl . [1989) has an inner torus photoionizing the geometrically thin 
disk that produces the disk line profile. However, the illuminated portion of the disk may have a 
complex structure. Studies of the variabil ity of double-peaked emitte rs have shown that there may 



be large-scale patterns in the emissivity ([Lewis and Eracleoust l2006h , suggesting that the surface 
emissi vity is not uniform . Proposed ideas for intrinsic asymmetries in ac cretion disks include hot 
spots ([Zheng et all Il99ll ) and regions of high density such as spiral arms ([Chakrabarti and Wiital . 



1994} that may result in a region of higher luminosity that is in the blue-shifted portion of the disk. 
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Such an asymmetry would produce the unusually strong peaks seen in some double-peaked Balmer 
emission lines. 



4. Conclusions 



In the case of SDSSJ0932+0318, the spectral features are unique and cannot be clearly classified 
as a typical double-peaked emitter. While maintaining that a disk profile is indeed present in the 
Balmer emission lines, we explored the scenarios of a chance projection, a binary SMBH, and a 
single accretion disk with a complex and asymmetric structure. The chance projection is unlikely, 
especially since there is only one set of narrow [OIII] lines. Since binary SMBHs are expected due to 
galaxy mergers, the model in which the strong blue peak in the disk is caused by a second SMBH 
is given serious consideration, especially since the velocity shift between the red and blue peaks 
in SDSSJ0932+0318 is similar to the shifts in the candidate binary SMBHs discovered previously 
and is on the high end for double peaked emitters. However, since much about the physics of 
accretion disks around SMBHs is not fully understood, we cannot discount the possibility of intrinsic 
asymmetries in the disk. 

Perhaps one of the best and easiest tests of the binary hypotheses is multiple observations, which 
may reveal orbital motion if it exists. From the basic binary parameters estimated in section 3.2, 
the hypothetical orbital period would be ~10 years. Such orbital motion may be apparent in 
the spectrum as a relative shifting of the two peaks over that time period, which would make it 
potentially observable. It has been pointed out that follow-up observations of t he three candidates 



previo usly mentioned (section 3.1) did n ot reveal a measur a ble velo city change (jBogdanovic et al 



2009al ). However, it was also noted by iBogdanovic et al] (|2009af) that this does not negate the 



binary SMBH hypothesis but may merely imply that the projected change in velocity will only be 
detectable over a period of several years or more. Even in the case where only a single SMBH 
is present, the strength of the blue H/3 and H7 emission peaks is notable and comparable to 
other unusual double-peaked emitters. Understanding the nature of the complex line emission in 
this object will aid in further discrimination between a single SMBH with a complex accretion 
disk and the actual case of a binary SMBH. Furthermore, the similarities of SDSSJ0932+0318 to 
other candidate binary SMBHs suggests that there may be a significant number of other objects 
with similar features that can be studied as binary SMBH candidates or unusual double-peaked 
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emitters. 
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